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Abstract 

We report properties of parametric electron pumping in the presence of a 



^- ■ superconducting lead. Due to a constructive interference between the direct 

reflection and the multiple Andreev reflection, the pumped current is greatly 
enhanced. For both quantum point contacts and double barrier structures 



at resonance, we obtain exact solutions in the weak pumping regime showing 
that Ip S = 4Ip , which should be compared with the result of conductance 



i-rt ■ Gns = 2Gtv. Numerical results are also provided for the strong pumping 

regime showing interesting Andreev assisted pumping behaviour. 
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Current can flow under zero bias when two system parameters of a nanostructure are 
varied in a cyclic fashion. The physics of this parametric electron pump has been analyzed 
by several authors [|1]|§. Recently, it has been realized experimentally by switkes et. al. || 
confirming many of the theoretical predictions. So far, investigations on parametric pumping 
are focused on normal nanostructures. It will be interesting to study a hybrid structure where 
a superconducting lead is present. In the presence of normal conductor-superconductor 
(NS) interface, an incoming electron-like excitation can be Andreev reflected as a hole-like 
excitation. This current doubling effect gives rise to the relation for conductance Gns — 
2Gn for hybrid quantum point contacts and quantum dots at resonance |4j]. For dirty NS 
contacts, de Jong and Beenakker || found that the shot noise at subgap voltages is also 
doubled with respect to its value in normal state which has been confirmed experimentally 
0. To further explore how does Andreev reflection modify quantum interference in the 
normal state, we have investigated the parametric pumping phenomenon in the presence of 
a superconducting lead. We find that due to quantum interference of the direct reflection 
and the multiple Andreev reflection, the pumped current is greatly enhanced. For quantum 
point contact and quantum dot at resonance, we obtained a relation for the pumped current, 
Ip = 4Ip in the weak pumping regime. Numerical results are presented in the strong 
pumping regime showing interesting Andreev assisted pumping behaviour which can be 
verified experimentally. 

We consider a parametric pump which consists of a double barrier tunneling structure 
attached to a normal left lead and a superconducting right lead. The double barrier structure 
is modeled by potential V(x) = Vi5(x + a) + V 2 (5(x — a) where V\ = V + V^sin(c<;£) and 
V2 = Vo + VpSm(ujt + (j)) and V p is the pumping amplitude. We further apply a gate voltage v g 
to control the energy level of the structure. The units are fixed by setting h = 1m = q = 1 in 
the following analysis |7| . At low frequencies, the adiabatic pumped current in the presence 
of superconducting leads is |1],|8| 



" " vrio T[ dV l dt dV 2 dt 



f = lMr[^ + ^] (1) 



where the quantity dNi/dV is the injectivity []9|JT0f] given, at zero temperature, by 

^ = ^MSledv, S ee - S^fy S he ] (2) 



where the first term is the injectivity of electron due to the external parameter PJTQ[|, i.e. 
the partial density of states (DOS) for an electron coming from left lead and exiting the 
system as an electron, and the second term is the injectivity of a hole, i.e. the DOS for a 
hole coming from left lead and exiting the system as an electron. 



For the hybrid nanostructure, the scattering matrix S ee and She are given by p|Jll 



S = Sn + 5i 2 (l — R1S22) R1S21 (3) 

where S is a 2 x 2 scattering matrix for NS structure with matrix element S^ (/i, v = e, h) 
and Sij(E) (i,j = 1,2) is a diagonal 2x2 scattering matrix for double barrier structure 
with matrix element Sij(E) and S*A—E). Ri = aa x is the 2x2 scattering matrix at 



NS interface with off diagonal matrix element a. Here a = (E — iv^A 2 — E 2 )/A with 
v — 1 when E > — A and v — — 1 when E < — A. In Eq.(^), the energy E is measured 
relative to the chemical potential \i of the superconducting lead. Eq.(BI) has clear physical 



meaning [II |. The first term is the direct reflection from the normal scattering structure 
and the second term can be expanded as S12R1S21 + S12R1S22R1S21 + ••• which is clearly the 
multiple Andreev reflection in the hybrid structure. From Eq.(|3|) we obtain the well known 
expressions for the scattering matrix S ee and She @ 

See(E) = S n (E) + a 2 S 12 (E)S*22(-E)M e S2i(E) (4) 

and 

She(E) = aS* 12 (~E)M e S2i(E) (5) 

with M e = [1 — a 2 822(E) S2 2 (~ ^)] ■ m the following, we first present the exact result 
for the pumped current in the weak pumping regime and then study the general situation 
numerically. In the weak pumping regime, Eq.fjj]) can be expanded to the lowest order in 
V p , 



NS u sin <pV t 



jn, = -_jL Im[dvi si e d V2 s ee - d Vl S* he dv 2 S he ] (6) 
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as compared with the expression for the normal structure jl|JT2 



AI u sin <t>V 2 
jN = T-JLlm^S'ndvzSn + d Vl S* 2l d V2 S 21 ] (7) 
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where we set V p = in S Ufl and Sy after the partial derivatives. We further assume that 
the Fermi energy is close to the chemical potential of superconducting lead, so E ~ and 
a ~ — i. Under this condition, She and hence dy She are pure imaginary numbers for general 
Vj(t). As a result, She does not contribute to Ip S in Eq.(^) as long as E = fTB| . The 
only contribution comes from S'ee which is the superposition of the direct reflection and the 
multiple Andreev reflection. We will consider two cases: (a) a quantum point contact, e.g. 
V = 0; and (b) the double barrier quantum dot at resonance. For both cases, Sn = in 
the absence of pumping voltage. Therefore, from Eqs.@ we have 

dv 1/2 S e e = d Vl/2 Su - S^dy^S*! (8) 



where we have used the fact that dy 1 S 2 2 = dy 2 Su. Using Fisher-Lee relation JT4 



Sa/3 = —S a /3 + i^/v^ypG^p and the Dyson equation d Vj G r a/3 = G r aj G r jf3 [0, we have 
dv 1 Sn = JvGijGy = —i/(2k) and dv 2 S\\ = %vG\ 2 G r 2X = —i/(2k)Sf 2 with the velocity 
v = 2k. So for both cases (a) and (b) we have constructive interference between direct 
reflection and the multiple Andreev reflection dy.S ee = 2dy .Sn which gives a pumped cur- 
rent Ip S = 41 p with 1^ = ±u)Vp sm(4ka)/(47rk 2 ), where the plus sign is for quantum 
point contact since S\ 2 = e 2lha , and the minus sign corresponds to resonant tunneling since 

C „—2ika 

<->12 — e 

In the general situation, the pumped current can be calculated numerically using Eq. ([!]) 
||16|| . Since the pumped current is proportional to u, we set oj = 1 for convenience. In the left 



inset of Fig. 1, we plot the ratio Ip S /Ip as a function of pumping strength V p /Vq for Vq = 79.2 
at the resonant point. We see that as the pumping amplitude increases, the constructive 
interference effect is suppressed. At small pumping strength, the ratio is about four which 



agrees with our theoretical analysis. At larger pumping strength, this ratio decreases to 
the value below two. Similar behaviour is seen for the quantum point contact. Fig.l shows 
the pumped current as a function of Fermi energy Ep for different pumping amplitudes. 
For Vo = 79.2, we have chosen v g = —9.39 so that one resonant level in the quantum dot is 
aligned with the chemical potential /x s of the superconducting lead in the absence of pumping 
voltage V p . For comparison, we also plot the Andreev reflection coefficient when V p = 0. 
Several observations are in order: (1). the pumped current is peaked near the resonant 
level showing clearly a resonant behaviour. This is because the pumped current (Eq.(|l|) 
is proportional to DOS of the system which reaches its maximum near the resonance. (2) 
As the pumping amplitude V p increases, the pumped current increases. (3). The pumped 
current has two asymmetric peaks. To understand this, we plot in the right inset of Fig.l 
the Andreev reflection coefficient versus Ep at different moments in one pumping cycle. 
From this inset, we observe that the Andreev reflection coefficient gives one or two peaks 
depending on the configuration of the system. This behaviour can be understood from 
the Breit-Wigner form of the resonant Anfreev reflection T& through a single level E = 
(measured relative to /x s ) [|T7|]: 

r 2 r 2 

Ta = L R (91 

A 4(E 2 + TA r /4y + T 2 L T 2 R K) 

where Ar = Tp — Tp and r = Tp + Tp. We see that Ta shows two peaks when Tp < Tp 
and just one peak otherwise. Note that in one pumping cycle half of the configurations 
corresponds to T L < T R , therefore two pumped current peaks show up in Fig.l because 
from Eq.(|l]) the pumped current is obtained through integral over all the configurations in 
one pumping cycle. Finally, the reason that two peaks are asymmetric is mainly due to 
the energy dependence of the self energy. If /i s is right in the middle of two resonant levels 
(Ei and E2), i.e., \i = (E\ + E2)/2, then electron coming from normal lead with incident 
energy E\ tunnels into the structure through the resonant level E\ and Andreev reflected as 
a hole back to the quantum dot through the resonant level E2 with a Copper pair created 
in the superconducting lead. We now examine the Andreev assisted pumping through two 



levels. For Vq = 79.2 and v g = —23.48, there are two resonant levels inside the subgap 
(fi s = 0) at Ei = 14.09 and E 2 = —14.09. Hence, strong Andreev reflections can occur 
near Ep = E\. Fig. 2 shows the two level pumped current versus Fermi energy (solid line). 
For comparison we also plot the corresponding Andreev reflection coefficient versus Fermi 
energy when the V p is switched off (see inset of Fig.2). Similar to Fig.l, the pumped current 
also shows strong resonant behaviour with smaller amplitude (compare Fig.l dot-dashed 
line). We found two peaks of pumped current around Ei, one is near the resonant energy 
and the other one is shifted to a smaller energy with larger current. Although it is similar to 
Fig.l but has different origin. This is due to the fact that when the pumping gate is turned 
on, the barrier heights and hence the resonant level change with time. To confirm this, we 
also plot the Andreev reflection coefficients T4 at several instants of one pumping cycle in 
the same figure. The peaks of T4 shift around the energy level at different pumping time 
and give the behaviour of pumped current. We have also calculated the pumped current for 
other system parameters and confirmed that the behaviour of pumped current shown here 
is generic. 

In summary, in the presence of superconducting lead, the pumped current is greatly en- 
hanced due to the quantum interference of direct reflection and multiple Andreev reflection. 
In the weak pumping regime, we have the relation Ip S = 41^ for both quantum point con- 
tact and the resonant tunneling structure. Interesting Andreev assisted pumping behaviours 
are revealed as well. 
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FIGURES 
FIG. 1. The pumped current versus Ep at different pumping amplitudes: V p = 0.03Vo, 

0.05Vo- The left inset: It; /It, versus V p /Vq at resonant point. The right inset: Andreev reflection 

coefficient Ta as a function of Ep with V p = 0.051/0 a t different pumping time: t = ir/4 (dot dashed 

line), t = tt/2 (dotted line), t = 37r/4 (solid line), and t = 7ir/4 (dashed line). System parameters 

in Fig.l and the insets: V = 79.2, v g = -9.39, <p = tt/2, and A = 20. 

FIG. 2. The two level pumped current I p as a function of Fermi energy (solid line) with 
Vp = 0.05Vo and v g = —23.481. We also plot Ta at several pumping time: t = tt/A (dotted line), 
t = tt/2 (dashed line), t = 37r/4 (dot-dashed line), and t = tt (long dashed line). The pumped 
current has been offset by 0.3 for illustrating purpose. Inset: Ta versus Fermi energy when V p = 0. 
Other parameters are the same as Fig.l. 
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